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ABSTRACT 



Context. Deep multicolour surveys are the main tool to explore the formation and evolution of the very faint galaxies 
which are beyond the spectroscopic limit with the present technology. The photometric properties of these faint galaxies 
are usually compared with current renditions of semianalytical models to provide constraints on the detailed treatment 
of the fundamental physical processes involved in galaxy formation and evolution, namely the mass assembly and the 
star formation. 

Aims. Galaxy counts over large sky areas in the 360nm near-UV band are particularly important because they are 
difficult to obtain given the low efficiency of near-UV instrumentation, even at 8m class telescopes. Observing in the 
near-UV bands can provide a first guess on the distribution of star formation activity in distant (up to z ~ 3) galaxies. 
A relatively large instrumental field of view helps in minimizing the biases due to the cosmic variance. 
Methods. We have obtained deep images in the 360nm U band provided by the blue channel of the Large Binocular 
Camera at the prime focus of the Large Binocular Telescope. We have derived over an area of ~ 0.4 sq. deg. the galaxy 
number counts down to [/ = 27 in the Vega system (corresponding to U=27.86 in the AB system) at a completeness 
level of 30% reaching the faintest current limit for this wavelength and sky area. 

Results. The shape of the galaxy number counts in the U band can be described by a double power-law, the bright side 
being consistent with the shape of shallower surveys of comparable or greater areas. The slope bends over significantly 
at [/ > 23.5 ensuring the convergence of the contribution by star forming galaxies to the Extragalactic Background 
Light in the near-UV band to a value which is more than 70% of the most recent upper limits derived for this band. 
We have jointly compared our near-UV and K band counts collected from the literature with few selected hierarchical 
CDM models emphasizing specific critical issues in the physical description of the galaxy formation and evolution. 

Key words. Surveys - Techniques: image processing - Galaxies: photometry - Galaxies: statistics 



1. Introduction 

Wide and deep multicolour surveys are useful tools to 
investigate in detail the processes of galaxy formation 
and evolution, especially beyond the spectroscopic ca- 
pabilities of current instrumentation. One of the main 
aims of the deep multicolour surveys is to provide a 
clear picture of the processes involved in the mass as- 
sembly and star formation of galaxies across the cosmic 
time. Galaxy observables, e.g. luminosity and mass 
functions, two point correlation function, are typically 
compared with current renditions of semi-analytical 
models in hierarchical CDM scenarios, in order to derive 
the key ingredients related to the physics of galaxy 
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The observed properties of the galaxy population, how- 
ever, are mainly affected by the lim ited statistics. Deep pen- 
cil beam su rveys, such as GOODS | Giavalisco et al. (2004) ) 
or HUDF ( Beckwith et al. (2003) i, have been carried out 
on relatively small sky areas and are subject to the cosmic 
variance effect. Larger surveys like COSMOS which extends 
over a 2 deg^ area are indeed shallower and limit the knowl- 
edge of the faint galaxy population at intermediate and high 
redshifts. 

To probe the statistical properties of the faint galaxy 
population reducing the biases due to the cosmic variance. 
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a major effort should be performed over large areas with 
efficient multicolour imagers at 8m class telescopes. This is 
especially true in the near-UV band (hereafter UV) where 
instrumentation is in general less efficient but where it is 
possible to extract information on the star formation activ- 
ity and dust absorption present in distant galaxies. 

In this context we are exploiting the unique 
power of the Large Binocular Camera (LBC) 
installed at the prime focus of the Large 
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allongo et al. 
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magnitude limits in the U band (A ~ 3600A) over areas of 
several hundreds of sq. arcmin. 

Long LBC observations in the UV band are of com- 
parable depth to those of the HDFs (although, obviously, 
the image quality will be poorer). Even a single pointing 
with LBC produces a field two orders of magnitude larger 
than that of the combined HDF-N and HDF-S. This is very 
important because the transverse extent of the HDFs cor- 
responds to about 1 Mpc at z 0.5 — 2 where Dark Matter 
clustering is still important. 

The goal of this paper is to provide UV (360nm) galaxy 
number counts down to the faintest magnitude limits avail- 
able from ground based observations. The comparison of 
normalization and shape of the observed counts with that 
predicted by theoretical hierarchical models can help to en- 
lighten critical issues in the description of galaxy forma- 
tion and evolution like dust extinction and the formation 
of dwarf galaxies. 

Throughout the paper we adopt the Vega magnitude 
system {Uvega = Uab — 0.86) and we refer to differential 
number counts simply as "number counts" , unless otherwise 
stated. 

2. The Data 

The deep UV observations described here have been 
carried out with the Large Binocular Camera (LBC, 



Giallongo et al. (2008)). LBC is a double imager installed 
at the prime foci stations of the 8.4m tele scopes LBT 
(Large Binocular Telescope, [Hill et aL (2008)[ ). Each LBT 
telescope unit are equipped with similar prime focus cam- 
eras. The blue channel (LBC-Blue) is optimized for imag- 
ing in the UV-B bands and the red channel (LBC-Red) for 
imaging in the VRIZY bands. The unvignetted FoV of each 
camera is 27 arcminutes in diameter, and the detector area 
is equivalent to a 23 x 23 arcmin^ field, covered by four 4K 
by 2K chips of pixel scale 0.225 arcsec. Because the mirrors 
of both channels are mounted on the same pointing sys- 
tem, a given target can be observed simultaneously over a 
wide wavelength range, improving the operation efhciency 
Extensive description of the twin LBC instrument can be 
found in Giallongo et al. (2008) Ragazzoni et al. (2006)[ 



Pedichini et al. (2003)[ JSpeziah et al (2008)[ 

We have used a deep U-BESSEL image of 3 hours, 
acquired in normal seeing condition (FWHM=1 arcsec) 
during the commissioning of the LBC-Blue camera, to 
derived faint UV galaxy-counts in a 478.2 arcmin^ 
sky area till U (Vega)= 26.5, in the Q0933+28 
field (Steidel et al. (2003)). The data have been re- 
duced using the LBC pipeline described in detail in 



age. The flux calibration of the U-BESSEL image has been 
derived through observations of photometric standards 
from the fields SA98 and SA1 13 (iLandolt et al. (1992) ) and 
the photometric fiel ds of iGaladi-e nriqucz ct al. (2000), as 



described in detail in Giallongo et al. (2008). The precision 
of the zero point calibration in the U-BESSEL filter is typ- 
ically of the order of 0.03 mag at 68% confidence level. A 
correction to the photometric zero point of 0.11 mag due to 
the Galactic extinction has been applied to the final coadd 
U-BESSEL image. 

The Q0933+28 field was also imaged in the SDT-UspeeQ 
filter of LBC for an additional hour in the first quarter of 
2007, under normal seeing conditions (1.1 arcsec). These 
images are reduced and coadded in the same way of the U- 
BESSEL ones, except for the photometric calibration pro- 
cedure, that is carried out using a spectrophotomctric stan- 
dard star of Oke (1990) The precision of the zero point 
calibration in the SDT-Uspec filter is typically of the order 
of 0.05 mag at 68% c.l. A Galactic extinction correction of 
0.13 has been used for this image. 

To push deeper the galaxy number counts in the UV 
band, we sum up the image obtained in the U-BESSEL fil- 
ter (with exposure time of 3 hours) with this new one, after 
rescaling the two images to the same zero point. We have 
verified that the effective wavelengths of these two filters are 
the same, the only difference being the higher transmission 
efficiency (1.5 times, after integrating its efficiency curve 
from A=3000 to 4000 A) of the SDT-Uspec filter compared 
to the U-BESSEL one. The colour term between these two 
filters is 0.01, thus we neglect it when summing up the two 
images. The resulting image goes '^0.5 mag deeper than the 
original 3 hours with U-BESSEL filter, given the higher im- 
age quality of the new SDT-Uspec image due to the general 
improvement of the telescope-instrument system, in partic- 
ular the reduction of scattered light from the telescope and 
dome environment after the first run of the LBC-Blue com- 
missioning. We use this final coadded image to improve the 
magnitude limit in the UV band and extend the number 
counts in this band to faint fluxes. 

To decrease the effects of cosmic variance in the num- 
ber counts at J7 ~ 20 we have used 3 additional LBC 
fields with shallower magnitude limits U < 25 but much 
larger area (892 arcmin?) in the Subaru XMM Deep Survey 



(SXDS, Sckiguchi et al. (2004) ) region. These images have 
been reduced and calibrated as described above for the 
Q0933-h28 field. The FWHM of these images (SXDSl, 
SXDS2, SXDS3) is higher (1.2, 1.25, 1.4 arcsec) than the 
one in the Q0933-|~28 field and the exposure time per LBC 
pointing is 1.0, 1.5, and 1.5 hours, respectively, in the U- 
BESSEL filter. These three LBC images, combined with 
the deep point in the Q0933-I-28 field are then used to de- 
rive wide and deep galaxy number counts in the U band, 
from U = 19.5 to 25.0 over a FoV of 1370 arcmin? (0.38 sq. 
deg.) and to ?7 = 27.0 for a sub-sample of 480 arcmir? . 

3. The number counts in the U band 

3.1. Deep U band galaxy number counts 

We com puted galaxy number co unts using the SExtractor 
package ( Bertin fc Arnouts (1996)| ). For objects with area 



Giallongo et al. (2008) applying the standard debias, flat- 
fielding and stacking procedures to derive the coadded im- 



this is an interference filter with a wavelength range similar 
to U-BESSEL filter, but the peak transmission is ^ 30% more ef- 



ficient than that, as described in Fig. 2 of Giallongo et al. (2008) 
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greater than that corresponding to a circular aperture of 
radius equal to the FWHM, we used the "best" photom- 
etry (Kron magnitude or corrected isophotal magnitude if 
the galaxy is severely blended with surrounding objects) 
provided by SExtractor. For smaller sources we computed 
magnitudes in circular apertures with diameter equal to 2 
times the FWHM, and correct them with an aperture cor- 
rection we derived using relatively bright stars in the field. 
This allows us to avoid the well known underestimate of 
the flux of faint galaxies provided by the isophotal method. 
To isolate the few stars from the numerous faint galaxies in 
this field, we relied on the class_star classifier provided by 
SExtractor. It is known that the morphological star/galaxy 
classifier of SExtractor is not reliable for faint objects, but 
the contamination from stars at faint U band fluxes and 
at high galac tic latitude is very limited (1%), according to 
the model of Bahcall & Soneira (1980) At brighter mag- 
nitudes, U ~ 20, the contamination from stars cannot be 
neglected, but at high S/N ratio the morphological clas- 
siflcator of SExtractor is robust. Moreover, the agreement 
between the LBC number counts at J7 < 22 with those of 
SDSS and other large area surveys (e.g., VVDS) indicates 
that the contamination from galactic stars are reduced also 
for bright U band magnitudes. 

We search for an optimal configuration of the detec- 
tion parameters of SExtractor in order to maximize the 
completeness at faint magnitude limits and to reduce the 
number of spurious sources. The two parameters regulating 
the depth, completeness and reliability of photometric cat- 
alog are the threshold (relative to the RMS of the image) 
adopted for detection of objects and the minimum area of 
connected pixels above this thres hold. We used the nega- 
tive image technique, described in Dickinson et al. (2004)[ 



Yan fc Windhorst (2004) 
timate of the reliability o: 



Bouwens~al. (2007) 



as an es- 
the catalog. Using SExtractor , 



we produce the "-OBJECTS" image, i.e. the image with 
the detected objects subtracted, then we compute the neg- 
ative of this image and run SExtractor with the same de- 
tection parameters used for the positive image. The sources 
detected on the negative image give an estimate of the con- 
tamination from spurious sources, provided that the noise 
statistic is symmetric (towards positive and negative pixel 
values) in the image. 

After several tests with different parameter set, the 
best effort combination are thresh=0.7, area=8 pixels for 
Q0933-h28 field and thresh=0.8, area=9 for SXDS fields 
(where the seeing is a bit worse) , which minimize the num- 
ber of spurious sources (measured in the negative image) 
while maintaining the number counts and completeness 
of real sources relatively high at faint magnitude limits 
{U ^ 26.5). In the Q0933-I-28 field the contamination from 
spurious sources is ~ 2.2% at U = 26.5, and remains less 
than 2.5% till U = 27.0 (corresponding to U(AB)=27.9). 

The results on the galaxy number count analysis can be 
affected by blending of galaxies when deep fields are inves- 
tigated. In SExtractor the two parameters affecting the de- 
blending of sources are the DEBLEND_NTHRESH and 
DEBLEND_MINCONT parameters. We produce differ- 
ent catalogs using various combination of these two pa- 
rameters for each LBC image, and we find that the raw 
number counts are not sensitive to these parameters, since 
the number of galaxies per magnitude bin and per square 
degree varies of ALogN = 0.02 at [/ = 24 - 25 and 
ALogN = 0.04 at U = 27, our faintest bin in the galaxy 



number counts. We adopt DEBLEND_NTHRESH=32 
and DEBLEND_MINCONT=0.002 for aU the fields de- 
scribed here, based on visual inspection of the reliability of 
deblended sources on the Q0933-I-28 field. 

Resulting raw counts are shown in FigU] where a clear 
decrease is apparent for U{Vega) > 26.4. The typical pho- 
tometric error at f7 ~ 27 is = 0.4 magnitude, correspond- 
ing to a total integrated number counts of ~100 galax- 
ies per arcmin? . An estimate of the completeness level 
should be performed in order to evaluate the amount of 
correction to the raw counts at the faint limits. This has 
been evaluated including in the real image 1000 simulated 
galaxies per 0.25 magnitude bin in the magnitude interval 
U{Vega) = 24 — 28 using the standard "artdata" package 
in IRAF. We include disk galaxies of half light radius of 
0.2-0.4 arcsec, convolved with the PSF of stellar objects 
in the field. We then run SExtractor on this new image us- 
ing the same detection parameters described above, and we 
studied how effective SExtractor is in recovering the sim- 
ulated galaxies as a function of magnitude. The resulting 
sizes of the simulated galaxies are typical of real galaxies 
in the magnitude interval U{Vega) = 24 — 25, so the com- 
pleteness we have computed at J7 = 26 — 27 can be con- 
sidered a rob ust estimate. We have ve rified using different 
datasets, i.e. [Windhorst et al. (2002)] and GOODS-South 



Grazian et al. (2006)[ that, at B > 25, the half light radius 



of these galaxies is always between 0.2 and 0.4 arcsec. We 
choose not to use real sources to simulate the completeness 
of the images mainly for two reasons: the half light radii 
of galaxies decrease with the magnitude, and it is known 
that fainter galaxies are smaller. If we use bright {U < 24) 
high signal-to-noise objects, their size, once dimmed to 
?7 ~ 26 — 27, can be overestimated with respect to their 
actual size, enhancing artificially the completeness correc- 
tion at fainter magnitudes. If we use real faint sources as 
input in our simulations, the signal-to- noise ratio of these 
galaxies is very low, enhancing even in this case the com- 
pleteness correction. We have tried instead a robust and 
realistic assumption, that is to simulate galaxies with half 
light radius between 0.2 and 0.4 arcsec, which is the range 
observed at the faint magnitude l imits where morphologi - 
cal analysis is still reliable (see e.g. Conselice et al. (2008) ). 
The completeness at U < 25.5 does not depend on the half 
light radius of the simulated galaxies, while at fainter mag- 
nitudes it depends on the size of the galaxies: in particular, 
at U — 26.0 the scatter of the completeness for various half 
light radii between 0.2 and 0.4 arcsec is 3%, increasing to 
6% atU = 26.5 and reaching 16% at U ^ 27.0. The result- 
ing 50% completeness level is measured at U (Vega) = 26.6, 
while at U = 27.0 we have a formal completeness of 30%. 

The number counts corrected for incompleteness are 
shown again in FigUJ Given the wide magnitude interval 
from U{Vega) = 19.5 to U{Vega) = 27.0 available in the 
present survey, the shape of the counts can be derived from 
a single survey in a self-consistent way, possibly minimiz- 
ing offsets due to systematics in the photometric analy- 
sis of data from multiple surveys (zero point calibration, 
field to field variation, etc). A clear bending is apparent at 
U{Vega) > 23.5. To quantify the effect we fitted the shape 
of the counts in the above magnitude interval with a double 
power-law. The slope changes from 0.58±0.03 to 0.24±0.05 
for magnitudes fainter than Ubreak = 23.6. The uncertainty 
in the break magnitude is however large, ~ 0.5, since the 
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transition between the two regimes of the number counts is 
gradual. 

In FigU] we compare our galaxy number counts 
with those derived by shallow s urveys of large area 



(SDSS EDR, 
(GOYA 
Radovich et al 



Yasuda et al. (2001) ) 
by |Etiche-M oral et a. 

(2004)p and with deep "pencil beam sur 



or 
"(^006) 



of similar area 
VVDS-F2 by 



veys (Hawaii HDFN by Capak et al (2004)[ WHT, HDFN, 
and HDFS by |Metcalte et al. (2001)p.~rn p articular, the 
WHT galaxy counts (Metcalfe et al. (2001)1 are based on 
a 34h exposure time image reaching U {V ego) = 26.8 but 
at the much lower 3cr level in the photometric noise and 
in an area of ~50 arcmin^, while the GOYA survey at the 
INT telescope is complete at 50% level at U{Vega) = 24.8. 
These counts are shown together w ith the two pencil beam 
surveys in the Hubble Deep Fields (Metcalfe et al. (2001) ). 
The agreement with the GOYA survey (900 sq. arcmin.) 
is remarkable, and suggests that once big areas of the sky 
are investigated, the effects of cosmic variance are slightly 
reduced. The present UV counts obtained during the 
commissioning of LBC-Blue are thus a unique combination 
of deep imaging in the U band and of large sky area, 
with the result of a considerable reduction of the cosmic 
variance effects for U > 21. For brighter magnitude limits, 
we refer to larger area surveys, shallower than our survey, 
as shown in FiglTJ 

Table [T] summarizes the galaxy number counts, cor- 
rected for incompleteness, with their upper and lower 1 
(T confidence level uncertainties, assuming Poisson noise 
and cosmic variance effect. For the latter, we used 
the Cosmic Variance Calculator (vl.02) developed by 



Trenti fc Stiavelh (2008) using as input values a linear size 
of 22 arcmin (corresponding to our deeper area of 478.2 
arcmin?) and redshift from 0.0 to 3.0, with standard A- 
CDM cosmology. At U=27, for example, the computed cos- 
mic variance is 4.5%. 

Using the Q0933-I-28 and the other three fields in the 
SXDS area, we study the field to field variation of the num- 
ber counts in the U band. We find that the typical variation 
from one LBC field to another is 0.04 in LogN for U ^ 20, 
while it reduces gradually to 0.01 at U=24.5, well below 
the poissonian uncertainties described in Table [TJ This en- 
sures that the zero-point calibration of the images is robust 
and the area of this survey is sufficient to decrease the cos- 
mic variance effects below the statistical uncertainties of 
the galaxy number counts. We complement these internal 
checks on the galaxy number counts in U with an external 
consistency test, comparing our final number counts with 
shallower NC derived in different areas or deep pencil beam 
surveys in Fig[TJ We find that the survey-to-survey maxi- 
mal variations are of the order of 0.1 in LogN at all magni- 
tudes, with lower scatter for the wider surveys considered 
here (GOYA, SDSS-EDR, VVDS-F2). 

3.2. The UV extragalactic background light 

The slope of the galaxy number counts in the U band at 
U ~ 23.5 changes from 0.58 to 0.24, which implies that the 
contribution of galaxies to the integrated EBL in the UV 
has a maximum around this magnitude. The contribution 
of observed galaxies to the optical extragalactic background 
light (EBL) in the UV band can be computed directly by 
integrating the emitted flux multiplied by the differential 
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Fig. 1. Number counts of galaxies in the U-BESSEL band 
for the Q0933-I-28 and SXDS LBC fields. Black crosses rep- 
resent the raw galaxy number counts, while filled big cir- 
cles show the number counts corrected for incompleteness. 
Magnitudes are in the Vega system. We compare our counts 
with shallow surveys of similar or larger area (SDSS EDR, 
GOYA, VVDS F2), and with deeper pencil beam surveys 
(Hawaii HDFN, WHT, HDFN, HDFS). 

Table 1. LBC U galaxy number counts. 



U(Vega) 


LogN 


Max LogN 


Min LogN 


Completeness 


19.50 


2.409 


2.648 


1.837 


1.00 


19.75 


2.484 


2.707 


2.001 


1.00 


20.00 


2.517 


2.734 


2.066 


1.00 


20.25 


2.891 


3.043 


2.654 


1.00 


20.50 


2.991 


3.129 


2.787 


1.00 


20.75 


3.103 


3.227 


2.930 


1.00 


21.00 


3.239 


3.346 


3.095 


1.00 


21.25 


3.350 


3.446 


3.227 


1.00 


21.50 


3.503 


3.585 


3.401 


1.00 


21.75 


3.655 


3.725 


3.572 


1.00 


22.00 


3.782 


3.843 


3.711 


1.00 


22.25 


3.914 


3.989 


3.824 


1.00 


22.50 


4.100 


4.161 


4.028 


1.00 


22.75 


4.191 


4.247 


4.127 


1.00 


23.00 


4.311 


4.360 


4.256 


1.00 


23.25 


4.452 


4.494 


4.405 


1.00 


23.50 


4.561 


4.598 


4.521 


1.00 


23.75 


4.633 


4.667 


4.595 


1.00 


24.00 


4.702 


4.734 


4.668 


1.00 


24.25 


4.799 


4.827 


4.768 


0.99 


24.50 


4.865 


4.892 


4.837 


0.98 


24.75 


4.945 


4.970 


4.919 


0.96 


25.00 


5.016 


5.040 


4.992 


0.93 


25.25 


5.081 


5.103 


5.058 


0.90 


25.50 


5.135 


5.156 


5.113 


0.88 


25.75 


5.199 


5.219 


5.178 


0.83 


26.00 


5.268 


5.288 


5.248 


0.75 


26.25 


5.341 


5.361 


5.321 


0.66 


26.50 


5.431 


5.450 


5.411 


0.55 


26.75 


5.469 


5.510 


5.425 


0.43 


27.00 


5.495 


5.563 


5.419 


0.30 



LBC U galaxy number counts corrected for incompleteness; A'^ is the 
number of galaxies per deg^ and per magnitude bin, while the minimum 
and maximum counts are the 1 cr confidence level (Poisson noise and 
cosmic variance effect). 
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Fig. 2. The extragalactic background light per magni- 
tude bin function of U band magnitudes. We com- 
pare the results of this w ork with the EBL derived by 
Madau fc Pozzetti (2000)] in the HDFS field. The smaU er- 
ror bars of EJiJL a.t U < 18 are due to the small uncertainties 
on the galaxy number counts derived using SDSS data. 



number counts down to the completeness limit of the sur- 
vey. ^ ^ 

' ' we 



Following the work of [Madau fc Pozzetti (2000) 
compute the EBL in the D band, 7^ measured in 
erg cm~^ Hz~^ Sr~^, using the following method: 

I, = 10-"-4(^-^«+4«-6)iV(t/Ai3) , (1) 

where we use the following correction from Vega to AB: 

UaB = Uvega + 0.86. 

Fig[2] shows the EBL obtained by the LBC deep 
number counts and compare it with the results of 



Madau fc Pozzetti (2000) Our data agrees well with pre- 
vious estimates of the HBL in the U band, and indeed it 
allows to derive precisely the peak of the differential EBL, 
at [/ = 23.55 ± 0.25, and put strong constraints to the 
contributions of faint galaxies down to U = 27.0 to the 
integrated EBL in this band. At magnitudes [/ > 24 our 
es timate of the EBL is sign ificantly larger than the one 
of [Madau fc Poz zet ti (2000)[ Checking in detail FiglH it is 
clear that the HDFB is slightly underdense with respect 
to the HDFN or other deep pencil beam surveys, probably 
due to cosmic variance effects, and this is reflected on the 



smaller EBL of Madau fc Pozzetti (2000) compared to the 
LBC estimate. 

The integrated galaxy contribution to the EBL in the 
UV band results = 3.428 ± 0M8nW/m^/Sr at 

the effective wavelength of the LBC U-BESSEL filter, 
3590A, taking into account the magnitude range U=17- 
27. This estimate is 20% h igher than the one estimated by 



Madau fc Pozzetti (2000)[ {2.87toilnW/m^ / Sr) although 



large uncertainties. Moreover, our estimate reduces the un- 
certainties by an order of magnitudes, which is of funda- 
mental importance to put strong constraints to galaxy evo- 
lution models. Extrapolating the observed number counts 
in the U band down to flux=0 we derive an integrated EBL 
of = 3.727 ± OMAnW/m'^ / Sr, so in our survey to 

U=27 we are resolving 92% of EBL produced by galaxies. 

The integrated flux from resolved galaxies, however, 
should be considered as a lower limit to the total EBL in the 
universe, since our estimate could be affected by different 
systematics: 

— Incompleteness in the number counts at faint mag- 
nitudes due to the presence of very extended/diffuse 
sources that escape detection in our LBC deep images 
and in the HDFs surveys. The adopted detection algo- 
rithms usually select galaxies down to a threshold in 
their surface brightness limits, and so are prone to se- 
lection effects against low surface brightness galaxies. 
We have computed the completeness of galaxies assum- 
ing a conservative half light radius of 0.3 arcsec, and 
we cannot exclude the presence of a population of dif- 
fuse and extended galaxies in the UV bands, even this 
hypothesis is quite unlikely, due to the well known rela- 
tion between obser ved optical magnitud es and half light 
radius of galaxies (Totani et al. (2001) ). 

— The surface brightness dimming of sources at high red- 
shifts. This is not a serious problem in the U band, since 
all the light in this wavelengths comes from ^ < 3 galax- 
ies (non U-dropout galaxies) and, as it is discussed in 
the next paragraph, the redshift distribution of galaxies 
contributing to the peak of the EBL are at z < 1. 



Independent estimates to the total E BL in the U band 

(2002)1 



given m 



are 

and are slightly 
galaxy number 



Bernstein et al. 
larger than 
counts {vl{v) 



our 



Bernstein (2007) 
estimate from 



14.4 ± 



Bernstein et al. (2002) [ and 



, idi^) = 

UA)nW/m''/Sr in [Bernstein (2007)[ after a 



tEe 

9nW/m^/Sr 
-- (21.6 ± 
different 



treatment for the contaminating foregrounds). 

W ith the present observa tions, we reject the hypothe- 
sis of Bernstein et al. (2002) that a huge contribution to 
the sky background in the U band could be derived by the 
overlapping wings of extended galaxies. With the present 
observations we reach a surface brightness limit oi U — 
2%. 2mag / arcsec^ at 3tT and thus exclude the presence of 
faint and extended wings in the galaxy UV light distribu- 
tions. 

The discrepancy between the resolved and total EBL 
could be explained by a population of ultra faint and nu- 
merous galaxies at [/ > 27, beyond the current limits of the 
present surveys and with a LogN-LogS slope steeper than 
0.5, or by an improp er subtraction o f t he foreground com - 
ponents, as s tated in |Mattila(2003)[ and [Bernstein (2007) 



Recently, Georganopoulos et al. (2008)[ proposed a new 



method to derive stringent limits to the diffuse light in the 
UV bands, using the inverse compton emission in 7-rays of 
the EBL by the high energy electrons in the radio lobes of 
Fornax A. This technique will allow to derive constraints to 



their value is still consistent with our measure due to their 



the TeV blazar emission 


(Stanev & Franceschini (1998) 


Aharonian et al. (2006) 




Albert et al. (2008)). At the 


present stage, Albert et a 

T/ \ ^ 1- ttA Tin 


. (2008) give an upper limit of 
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limit strengthens the hypothesis that the discrepancy 
between our estimate of the EBL in U band and that of 



Bernstein (2007) is due to a foreground local component. 

Our derivation of the EBL by integrated galaxy counts 
is marginally consisten t with the lower limit derived by 
Kneiske fc Dole (2008)| of > 3.97nW/rn^ / Sr at A = 

3600A. 

4. Comparison with theoretical models 

The availability of galaxy number counts down to 
faint magnitude limits and over large sky areas can 
be used to test the predictions of different theoreti- 
cal models without being strongly affected by the cos- 
mic variance. Among the various models (numerical 
and semi-analytical) developed in the framework of the 
standard hierarchical CDM scenario, we have selected 



the mode ls developed by Menci et al. (2006) (hereafter 
M06), the [Kitzbichler fc White (2007) model based on the 
Millenniu m Simulation (hereafter K07) and the model de- 
veloped by [Monaco et al. (2007), (hereafter MORGANA). 
All these model simulate the formation and evolution of 
galaxies, starting from a statistical description of the evo- 
lution of the DM halo population (merger trees), and us- 
ing a set of approximated, though physically motivated, 
"recipes" to treat the physical processes (gas cooling, star 
formation, AGN and SF feedback, stellar population syn- 
thesis) acting on the barionic component. An appropri- 
ate treatment of dust attenuation is also crucial in com- 
paring model predictions to the UV observations, given 
the efficient scattering of radiation at these wavelengths 
by dust grains, as shown by GALEX number counts of 
Xu et al. (2005)] 

The comparison of the observed counts with that pre- 
dicted by the three selected models is used to enlighten 
different critical issues concerning the physical description 
of the galaxy formation and evolution. 

In particular, the M06 model is characterized by a spe- 
cific implementation of the AGN feedback on the star for- 
mation activity in high redshift galaxies. The description 
is based on the expanding blast waves as a mechanism to 
propagate outwards the AGN energy injected into the in- 
terstellar medium at the center of galaxies (see Menci et 
al. 2008); such a feedback is only active during the active 
AGN phase ( "QSO mode" ) , and it is effective in suppress- 
ing the star formation in massive galaxies already by z ~ 2, 
thus yielding a fraction of massive and extremely-red ob- 
jects in approximate agreement with observations (Menci 
et al. 2006). The M06 model uses three different prescrip- 
tions for the dust absorption, namely the Small Magellanic 
Cloud model (SMC) the Milky Way (MW) law and the 
Calzetti et al. (2000)| extinction curve (COO). In FigE] the 
two short-dashed curves enclose the minimum and maxi- 
mum number counts predicted by the M06 model taking 
into account the three different extinction curves (SMC, 
MW, and COO). 

The model of [Kitzbichler fc White (2007)1 is ba sed on 
the Millennium Simulation (Springel et al. (2005)1, on a 
concordance Lambda CDM cosmology. In the model an im- 
portant critical issue is represented by the particle mass res- 
olution of the code, which affects the physical and statisti- 
cal properties of faint low-luminosity galaxies. The adopted 
CDM mass resolution is sufficient to resolve the halos host- 
ing galaxies as faint as 0.1 L* with at least 100 particles 



of 8.6 ■ IO^H-^Mq. The model also adopts a different AGN 
feedback called " radio mode" in which cooling is suppressed 
by the continuous accretion of hot gas onto Supermassive 
Black Holes at the center of groups or cluster of galax- 
ies. The only difference of K07 from Croton et al. (2006) 
and De Lucia & Blaizot (2007) is in the dust model used in 
their simulations. I'or the local galaxies, they adopt a sim- 
ple relationship between face-on optical depth and intrinsic 
luminosity r oc ro(L/L*)^ with (3 = 0.5. At high redshifts 
they take into account the different dust and gas contents, 
the varying metallicities, and the shorter rest frame emitted 
wavelengths. In particular their average extinction increases 
strongly at high red shifts due to the smaller disc s izes of the 
galaxies. We refer to Kitzbichler fc White (2007) for details 
about their dust model. 

Finally, the MORGANA model is characterized by a 
different treatment of the processes of gas cooling and infall 
(following Viola et al., 2008), star formation and feedback 
(using the multi-phase model of Monaco 2004) . Black hole 
accretion is described in detail in Fontanot et al., (2006): 
the accretion rate in Eddington units determines the nature 
of feedback from the AGN. Synthetic SEDs for model galax- 
ies are obtained using the GRASIL code (Silva et al. 1998), 
which explicitly solves the equation of radiative trans- 
fer in a dusty medium. Fontanot et al. (2007) show that 
MORGANA is able to reproduce both sub-mm (850/im) 
number counts and the redshift distribution observed in 
if-limited samples, with conservative choice for the stellar 
IMF. Galaxy SEDs, magnitudes, and colours for a variety 
of passbands are obtained using the GRASIL spectropho- 
tometric code (Silva et al. 1998), which explicitely solves 
the equations for the radiative transfer in a dusty medium 
taking into account the composite geometry (bulge-|-disc) 
of each model galaxy. The dust properties (dust-to-gas 
mass fractions, composition, and size distribution) are kept 
fixed to values providing a good agree ment with local ob- 
servations (see Fontanot et al. (2007) for more details on 
the couphng between MORGANA and GRASIL). In this 
work we take advantage of the recent update of the model, 
adapted for WMAP3 cosmology and Chabrier IMF and pre- 
sented in Lo Faro et al. (2009). 

It is worth noticing that photoelectric absorption by the 
interstellar medium is redshifted, for high redshift galaxies, 
into the observed UV band. For this reason only galaxies 
with z < 3 contribute to the UV number counts. The com- 
mon characteristics of the models is that only galaxies in 
the 1.5 < z < 2.5 redshift interv als give the main co ntribu- 
tion to the counts at J7 > 26 (see Barro et al. (2009)). This 



implies that the shape of the faint counts is provided by the 
average faint-end shape of the galaxy luminosity function 
in the same redshift interval. 

The comparison of the observed galaxy number counts 
in the U (360 nm) band with the predictions of theoretical 
models shown in FigO indicates that: 

— the M06 and K07 models reasonably reproduce the UV 
number counts from U = 17 to U = 27, while the 
MORGANA model shows the largest deviation from 
the observed data. 

— The three models show different shapes for the UV 
galaxy number counts at faint fluxes, implying discrep- 
ancies which increase for increasing magnitudes. This 
can be evaluated in more detail divid ing the counts in 
redshift bins (see Barro et al. (2009)). 
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Fig. 3. Number counts of galaxies in the U-BESSEL band 
for the Q0933+28 LBC field, complemented at brighter 
magnitudes {U < 20) by the number counts of VVDS-F2 
and SDSS-EDR, whose error bars are always smaller than 
the size of the points. We compare our counts with theo- 
retical models (M06, K07, MORGANA). 
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Fig. 4. Number counts of galaxies in the K band for differ- 
ent deep surveys in the literature. The model predictions 
come from the same models described in the previous fig- 
ure. 



— With the current observations of the UV number counts 
down to ?7 = 27 both M06 and K07 models show a 
similar behaviour, although the K07 model appears to 
bend over for [/ > 27 in contrast with the steeper slope 
of the M06 model. 

The discrepancies between predicted and observed faint 
UV galaxy number counts could be due in principle to 
several effects, e.g. the evolution in number density, star 
formation activity or dust extinction of the faint galaxy 
population. Among these, an appropriate treatment of the 
dust extinction plays an important role, in this comparison 
based on the U band, since it can affect both the normaliza- 
tion and shape of the observed UV counts. The effects due 
to the different extinction laws (Calzetti 2000 and Small 
Magellanic Cloud) adopted in the M06 model for example, 
increase at fainter magnitudes since we are observing galax- 
ies at high {z < 3) redshifts. In this respect, the bending 
observed in the UV counts could be indicative of a gradual 
redshift evolution of the dust properties. 

Another physical quantity relevant for the correct re- 
production of the faint counts is the amount of star for- 
mation activity in faint low-mass galaxies. This effect is 
particularly important in the MORGANA model where it 
is responsible for the flat shape of the predicted counts. 
Fontanot et al. (2009) show that the underprediction of UV 
number counts of tSnt galaxies can be explained by a rapid 
decline in their star formation activity and consequently in 
their associated UV emission at z < 2. These galaxies are 
predicted to be too passive and to host too old stellar pop- 
ulations at later times with respect to observations. The 
reduced SFRs can easily explain the flat shape of the pre- 
dicted counts in MORGANA. 



As an attempt to disentangle dust/SFR effects from 
number density evolution, we have extended the model 
comparison to the NIR K band number counts which 
are much less sensitive to dust absorption and short 
episodes of star formation. Indeed, the observed flux 
in the K band is more related to the star formation 
history in the galaxy quantified by its assembled stel- 
lar mass (see e.g. Fontana et al. 2006). Fig[4] shows 
the comparison of the observed number counts in the 
K band collected in the literature and provided by 
different instrument s / telescop es and surveys (SSDF: 
'Minowa et al. (2005)[ H DFS: ILabbe et al. (2003)i KDS: 
Moustakas et al. (1997)| UPS: [Cirasuolo et al. (2008)[ 
WHTDF: Cristobal-Hornillos et al. (2003)[ HWDF: 



Huang et al. (1997), GOODS: ^azian et al. (2006)p with 
'the predictions of M06, K07, and MORGANA. The 
observed K band galaxy counts show a clear bending at 
K 2i 17 where the average slope changes from 0.69 to a 
flatter value 0.33. Here the bending is more clear compared 
to what observed in the UV since the break magnitude is 
brighter with respect to the faintest survey limits. 

In the K band the models tend to overestimate the num- 
ber counts both at the bright end if < 15 and at the faint 
magnitudes K > 20. The excess in the predicted K band 
number counts resembles that at the faint end of the the- 
oretical luminosit y fu nctions in the NIR (se e for example 



PoU et al. (2003)1 and [Fontanot et al. (2009) ). This can be 
indicative of an excess in the production of galaxies with 
low luminosity and high redshifts. If the NIR number counts 
are extrapolated at fainter limits, the amount of this excess 
differs among the models, being smaller in the K07 model 
at K > 26. The latter however is affected by a thresh- 
old in the mass resolution used by the numerical code of 
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the Millennium simulation adopted in the K07 model. This 
could provide an artificial removal of faint galaxies. 

The simultaneous comparison of the UV and NIR K 
band number counts at faint fluxes seems to indicate that 
the underprediction of the modelled UV number counts 
could be due to dust extinction or reduced SF activity 
rather than to an intrinsic evolution in the galaxy number 
density. 

Finally we identify a third issue when comparing the 
bright end of the galaxy counts. The treatment of the feed- 
back on SF processes due to AGN activity can play an 
important role in this respect. The M06 model for exam- 
ple shows the effect of a different treatment of AGN feed- 
back based on the so called QSO mode at variance with 
the Radio mode feedback used by the K07 model, while 
MORGANA is characterized both by Radio and by QSO 
modes. The different treatment of AGN feedback is plau- 
sibly the responsible for the differences in the predictions 
of NIR galaxy number counts at K < 15. The excess of 
the UV number counts predicted by models is related to 
the quenching of the star formation activity. MORGANA 
tend to overpredict star formation in massive central galax- 
ies at low-z. This is related to the less efficient, or delayed, 
quenching of the cooling flows in massive halos via Radio- 
mode feedback. In fact, in this model gas accretion onto 
the central black hole is related to star formation activity 
in the spheroidal component: this implies that AGN heating 
switches on only after some cooled gas has already started 
forming stars in the host galaxy (see Kimm et al., 2008 for 
a complete discussion and a comparison of different Radio- 
mode implementations in semi-analytical models). 

As a last comment, all the three models give a contri- 
bution to the UV EBL which is broadly consistent with 
the upper limits described in section 3.2. Using the same 
method described there for the observed data, we derived 
an EBL of 0.71-1.18 nW/m'^/Sr for the M06 model in the 
UV band, while the MORGANA and K07 ones give 2.66 
and 3.26 nW/m? / Sr, respectively. 

5. Summary 

To summarize the main results of the paper: 

— We have derived in a relatively wide field of 0.4 deg^ 
the deepest counts in the 360nm UV band. This allowed 
to evaluate the shape of the galaxy number counts in a 
wide magnitude interval U= 19-27 with the advantage of 
mitigating the cosmic variance effects. The agreement 
with the number counts of shallower surveys confirms 
the low impact of systematic errors on the LBC galaxy 
statistics. 

— The shape of the counts in UV can be described by a 
double power-law with a steep slope 0.58 followed be- 
yond U ~ 23.5 by a flatter shape 0.24. Our counts are 
consistent at the bright end with surveys of comparable 
or greater areas. At the faint end our counts are more 
consistent with that found in the HDF-N. 

— The faint-end slope of the counts is below 0.4 and 
this ensures the convergence of the contribution by 
star forming galaxies to the EBL in the UV band. 
The total value in the UV band obtained extrapo- 
lating the slope of our counts to flux=0 is indeed 
3.727 ± 0.084nW/m'^/Sr. It is consistent with re- 
cent upper limits coming from TeV observations of 



Albert et al. (2008)[ < 5nW/m^/Sr, showing 

that the U V EBL is resolved at >74% level. 

— We have compared our counts in the UV and K bands 
with few selected hierarchical CDM models which are 
representative of specific critical issues in the physical 
description of the galaxy formation and evolution. 

— The mass resolution of numerical models is critical for 
reproducing the faint end of the UV galaxy number 
counts. 

— The discrepancies between predicted and observed UV 
galaxy number counts at faint magnitudes could be 
mainly due to the treatment of dust extinction and the 
star formation activity in low mass galaxies at z < 2. 

— The AGN feedback (Radio vs QSO mode) may affect 
galaxy counts at the bright end of the LogN-LogS in 
the K band. 

A correct physical description of the AGN feedback, 
dust properties and star formation activities in the mod- 
els is fundamental to ensure a reasonable agreement of the 
model predictions at the faint end of the galaxy counts. 

Adding colour information for galaxies with UV emis- 
sion as faint as f7 = 27 — 28 implies very deep observations 
in the red bands which are feasible with several hours of 
integration at 8m class telescopes. Very deep multicolour 
information on areas of the order of the square degree can 
help in extracting physical information on the star forma- 
tion history of the dwarf population at intermediate and 
high redshifts. 

Acknowledgements. Observations have been carried out using the 
Large Binocular Telescope at Mt. Graham, Arizona, under the 
Commissioning phase of the Large Binocular Blue Camera. The LBT 
is an international collaboration among institutions in the United 
States, Italy and Germany. LBT Corporation partners are: The 
University of Arizona on behalf of the Arizona university system; 
Istituto Nazionale di Astrofisica, Italy; LBT Beteiligungsgesellschaft, 
Germany, representing the Max-Planck Society, the Astrophysical 
Institute Potsdam, and Heidelberg University; The Ohio State 
University, and The Research Corporation, on behalf of The 
University of Notre Dame, University of Minnesota and University 
of Virginia. The Millennium Simulation databases used in this pa- 
per and the web application providing online access to them were 
constructed as part of the activities of the German Astrophysical 
Virtual Observatory. Some of the calculations were carried out on 
the PIA cluster of the Max-Planck-Institut fiir Astronomie at the 
Rechenzentrum Garching. We thank the anonymous referee for useful 
comments which helps in improving the quality of the present pa- 
per. AG warmly thanks Kalevi Mattila and Martin Raue for useful 
comments on the EBL limits. 



References 

Aharonian, F., et al. 2006, Nature, 440, 1018 



Albert, J. et al. 2008, Science, 320, 1752; arXiv:0807.2822 [astro-ph] 
Bahcall, J. N., & Soneira, R. M., 1980, ApJS, 44, 73 
Barro, G., Gallego, J., Perez-Gonzalez, P. G., et al. 2009, A&A, 494, 
63 

Beckwith, S. V. W., et al. 2003, AAS, 202, 1705 

Bernstein, R. A., Freedman, W. L., Madore, B. F., 2002, ApJ, 571, 
56 

Bernstein, R. A., 2007, ApJ, 666, 663 

Berlin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Bower, R. G., Benson, A. J., Malbon, R., et al. 2006, MNRAS, 370, 
645 

Bouwens, R. J., lUingworth, G. D., Franx, M., Ford, H. 2007, ApJ, 
670, 928 

Calzetti, D., Armus, L., Bohlin, R. C, et al. 2000, ApJ, 533, 682 
Capak, P., Cowie, L. L., Hu, E. M., et al. 2004, AJ, 127, 180 
Cirasuolo, M., McLure, R. J., Dunlop, J. S., et al. 2008, 
iarXiv:0804.3471 [astro-ph] 



Grazian et al.: Wide and deep near-UV galaxy counts with LBC 



9 



Conselice, C. J., Rajgor, S., Myers, R., 2008, MNRAS, 386, 909 
Cristobal-Hornillos, D., Balcells, M., Prieto, M., et al. 2003, ApJ, 595, 
71 

Croton, D. J., et al. 2006, MNRAS, 365, 11 

Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451 
De Lucia, G., Blaizot, J., 2007, MNRAS, 375, 2 
Dickinson, M., et al. 2004, ApJ, 600, L99 

Eliche-Moral, M. C., Balcells, M., Prieto, M., et al. 2006, ApJ, 639, 
644 

Fontana, A., Salimbeni, S., Grazian, A., et al. 2006, A&A, 459, 745 
Fontanot, F., Monaco, P., Cristiani, S., & Tozzi, P. 2006, MNRAS, 
373, 1173 

Fontanot, F., Monaco, P., Silva, L., Grazian, A., 2007, MNRAS, 382, 
903 

Fontanot, F., De Lucia, G., Monaco, P., Somerville, R. S., Santini, P. 

2009, arXiv:0901.1130 [astro-ph] 
Galadi-enriquez, D., TruUols, E., Jordi, C. 2000, A&AS, 146, 169 
Georganopoulos, M., Sambruna, R. M., Kazanas, D., et al. 2008, ApJ, 

686, L5 

Giallongo, E., Ragazzoni, R., Grazian, A., et al. 2008 A&A, 482, 349 
Giavalisco, M. and the GOODS Team 2004, ApJ, 600, L93 
Grazian, A., Fontana, A., De Santis, C., et al., 2006, A&A, 449, 951 
Hill, J. M., Green, R. F., Slagle, J. H., et al., 2008, Proc S.P.LE., 7012, 
pp.701203-15. 

Huang, J. S., Cowie,, L. L., Gardner, J. P., et al. 1997, ApJ, 476, 12 
Keres, D., Katz, N., Dave, R., Fardal, M., Weinberg, D. H., 2009, 

arXiv:0901.1880 [astro-ph] 
Kimm, T., Somerville, R. S., Yi, S. K., et al. 2009, MNRAS, 394, 1131 
Kitzbichler, M. G., & White, S. D. M., 2007, MNRAS, 376,2 
Kneiske, T. M. & Dole, H., 2008, High Energy Gamma-Ray 
Astronomy: 4th International Symposium, Proceedings of the 
conference held July 7-11, 2008 in Heidelberg (Germany) 
Labbe, I., Franx, M., Rudnick, G., et al. 2003, ApJ, 125, 1107 
Landolt, A. U. 1992, AJ, 104, 372 

Lo Faro, B., Monaco, P., Vanzella, E., et al. 2009 in prep. 
Madau, P. & Pozzetti, L. 2000, MNRAS, 312, L9 
Mattila, K. 2003, ApJ, 591, 119 

Menci, N., Fontana, A., Giallongo, E., Grazian, A., Salimbeni, S., 

2006, ApJ, 647, 753 
Metcalfe, N., Shanks, T., Campos, A., McCracken, H. J. & Fong, R. 

2001, MNRAS, 323, 795 
Minowa, Y., Kobayashi, N., Yoshi, Y., et al. 2005, ApJ, 629, 29 
Monaco, P. 2004, MNRAS, 352, 181 

Monaco, P., Fontanot, F., Taffoni, G., 2007, MNRAS, 375, 1189 
Moustakas, L. A., Davis, M., Graham, J. R., et al. 1997, ApJ, 475, 
445 

Nagamine, K., Ostriker, J. P., Fukugita, M., Gen, R. 2006, ApJ, 653, 
881 

Oke, J. B., 1990, AJ, 99, 1621 

Pedichini, F., et al. 2003, Proc. SPIE, 4841, 815 

Poh, F., Giallongo, E., Fontana, A., et al. 2003, ApJL, 593, 1 

Radovich, M., Arnaboldi, M., Ripepi, V. et al., 2004, A& A, 417, 51 

Ragazzoni, R., et al. 2006, Proc. SPIE, 6267, 33 

Sekiguchi, K.; et al. 2004, Multiwavelength Cosmology. Proceedings of 
the "Multiwavelength Cosmology" conference, held on Mykonos 
Island, Greece, 17-20 June, 2003. Ed. M. PUonis. 

Silva, L., Granato, G. L., Bressan, A., & Danese, L. 1998, ApJ, 509, 
103 

Somerville, R. S., Hopkins, P. F., Cox, T. J., et al. 2008, MNRAS, 
391, 481 

SpeziaU, R., Pedichini, F., Di Paola, A., et al. 2004, SPIE, 5492, 900 
SpeziaU, R., Di Paola, A., Giallongo, E., et al. 2008, SPIE, 7014, 158 
Springel, V., et al. 2005, Nature, 435, 629 
Stanev, T. & Franceschini, A., 1998, ApJ, 494, L159 
Steidel, C. C, Adelberger, K. L., Shapley, A. E., et al. 2003, ApJ, 
592, 728 

Trenti, M. & Stiavelli, M. 2008, ApJ, 676, 767 

Totani, T., Yoshii, Y. Iwamuro, F., et al. 2001, ApJ, 550, L137 

Viola, M., Monaco, P., Borgani, S., Murante, G., & Tornatore, L. 

2008, MNRAS, 383, 777 
Windhorst, R., Cohen, S., Jansen, R., et al. 2002, AAS, 201, 3207 
Xu, C. K., Donas, J., Arnouts, S., et al. 2005, ApJ, 619, Lll 
Yan, H. & Windhorst, R. A. 2004, ApJ, 612, L93 

Yasuda, N., Fukugita, M., Narayanan, V. K., et al., 2001, AJ, 122, 
1104 



